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Abstract 



A calculation of deuteron Compton scattering using non-relativistic pertur- 
\ bation theory is presented, with the primary motivation of investigating the 

I feasibility of determining the neutron polarizabilities from this type of exper- 

, , ^ iment. This calculation is expected to be valid for energies below 100 MeV. 

^ I Pion-exchange, relativistic, and recoil corrections are also included. The low- 

energy theorem for gauge invariance is shown to be satisfied. The relative 

><■ 

5— ( ' effects of the different terms and their effects on the determinations of the po- 

\ 

larizabilities are discussed at energies of 49, 69, and 95 MeV. The cross-section 
is dominated by the seagull, polarizability, and electromagnetic multipole in- 
teractions. Relativistic and pion-exchange terms are also important, while 
recoil corrections and multipoles of L=2 and greater are negligible. The cal- 
culation provides a reasonable description of the experimental data points at 
49 and 69 MeV. The polarizabilities are difficult to determine at these ener- 
gies. A more accurate determination of the polarizabilities may be possible 
at 95 MeV. 
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I. INTRODUCTION 



This paper is concerned with determining the polarizabihty of the neutron using Compton 
scattering from a deuteron. This is not a new idea and there have been previous calculations. 
However, advances on both the theoretical and experimental fronts invite the possibility of 
not only a more accurate calculation, but a more accurate measurement. The intent of this 
work is provide a complete potential model calculation. 

A deuteron Compton scattering experiment was performed in 1994 using 49 and 69 
MeV photons, and experiments are in progress at Saskatoon |^ and Lund This has 
spurred recent activity. Brief articles by Wilbois, Wilhelm, and Arenhovel, [Q and L'vov 
and Levchuk |^,|[ quote new theoretical results for Compton scattering. Calculations using 
effective field theory, based on the approach of Kaplan, Savage and Wise , have also been 
performed Another calculation [0, using an approach intermediate between ours and 
that of Ref . ||,^ , is nearly complete. 

Here is an outline of this paper. The remainder of this introduction reviews the concept 
of polarizabihty and examines the previous theoretical and experimental work. Descrip- 
tions of the methods used in the calculation are discussed in Sect. II. Checks which were 
performed to ensure the correctness of the work are also detailed there. The results of the 
numerical calculation are presented in Sect. HI. The possibility of determining the neutron 
polarizabilities from deuteron Compton scattering is detailed. The conclusion and summary 
are in Sect. IV. Comparisons between our and previous calculations will be made in the 
summary. Some of the details of calculations omitted in the main text are in the Appen- 
dices. Complete formulae for all of the relevant scattering amplitudes may be found in the 



1998 Ph. D. thesis 111 of one of the authors. This work is based on that thesis. 
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A. Background and Motivation 



The symbol a is traditionally used to represent the electric polarizability, and (5 is used for 
the magnetic polarizability. The polarizability is the constant of proportionality between an 
external field and the average dipole moment that it induces. This induced moment in turn 
changes the potential energy of a system by an amount IS.E = —^aE"^ — \l3B'^, known as the 
quadratic Stark effect. The nucleon polarizabilities are obtained using Hamiltonians for a 
dipole interacting with an external field: i^int = —d^E, or —fi^B. Second order perturbation 
theory leads to 



« = 22^ — p p , P = ^2^ — p p • (1) 



We make an order-of-magnitude estimate of a by letting A^' be an A^* state of mass 1440 
MeV and taking the square of the matrix element to be about |fni^. Then a ~ 13 x 10~^fm^. 
This is of the same order as the currently accepted value of a ~ 12 x 10~^fm^. Hereafter 
the unit of polarizability will be 10~^fm^ . 

The above expressions define the polarizability in terms of either an energy shift or 
an induced dipole moment. A third definition, which turns out to be the most practical, 
involves the scattering amplitude for elastic photon scattering, commonly known as Compton 
scattering. An initial photon, of laboratory four-momentum (cUj, fcj) and polarization vector 
exii where Aj = ±1, is absorbed by a nucleon, which emits a final photon of four-momentum 
{ujf, kf) and eA^. The angle between ki and kf is labeled 6, and the initial and final photon 
energies are related by Uf = oji/ (1 + — cos 6)). The scattering amplitude (to order cj^) 



IS given as 
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f{uj,uji) = (et -exj + {uji + u f)g{e,mN, kn] kf, el ki, ex,) + 

aujfUJi{el^ ■exd+ P^fUJi{kf ^ ^*Xf) ' x e^) + 

ujfUJih{e, mjv, k^] kf, e*x^,ku h,)- (2) 

The first term, the Thomson amplitude, is the dominant one at low energies. The internal 
structure becomes important only in second order. The terms with the angular dependence 



(^Aj ■ ^aJ and (kf x e^^,) • (fcj x e\.) are written separately, and all other combinations of the 
four photon unit vectors are lumped together in the function h. The coefficients of these 
two terms are defined as the generalized polarizabilities a and f5. 

The generalized polarizabilities a, (3 are not the same as the static polarizabilities a,(3. 
They are related by a = a + Aa, (3 = (3 + A(3, where Aa, Af3 arise from terms in the 
scattering amplitude of the same energy and angular dependence as those of a, (3. The first 



term in the explicit expression for Aa is [ll3| , p^ : Aa = + ^-j-^^^, where re stands 



3m 4m'^ 

for the charge radius, and accounts for about 40% of a^. This is only about a 5% effect in 
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the neutron. The first term in A/S is A/3 = — Additional corrections involve the 
nucleon form factors Fi(4m^), F2(4m^) and their derivatives M. 



Our primary interest is in determining polarizabilities through experiment, but it is useful 
to review the theoretical calculations. Baldin obtained an estimate 4.0 < a^ < 15.0, 
for the proton polarizabilities by including the first excited state ( A^' = N + tt). Gell-Mann 
and Goldberger |jT6| derived the once-subtracted dispersion relation for a forward scattering 



amplitude f{uj). The low energy limit of Eq. (|^), along with the dispersion relation and the 



the optical theorem, leads to the Baldin sum rule |15 



where (Ttot('^) is the total cross-section for photoabsorption. This relation provides a model- 
independent constraint on the possible values of a and /?. The generally accepted results 
from the sum rule are [p!7| , p!8| a^ + /^p = 14.3 ± 0.5, a„ + /?„ = 15.8 ± 0.5. However, a recent 
experiment [|19| predicts: a^ + /?p = 13.69 ± 0.14, a„ + /?„ = 14.44 ± 0.69. The values for the 



proton are consistent with the older ones, but the new neutron values are slightly lower. 
A dispersion relation for a — (3 can also be derived pO|j2T[] , but it contains a contribution 



that depends on the amplitudes for 77 ^ tttt, which are not well known. Nevertheless, this 



difference has been estimated to be ap — /3p ~ 3.2, a^ — /3„ ~ 3.9. Comparison with the 
sum rules above yield smaller values for the individual polarizabilities than the ap ~ 11 and 
a„ ~ 12 measured experimentally. 



Various quark models have also been used. A simple nonrelativistic model with a har- 
monic oscillator potential yields reasonable values for (3 but predicts > q;„, so more 
sophisticated models are needed. For example, Werner and Weise obtained a and (3 
using a valence quark core surrounded by a pion cloud to be ~ 7 — 9, /^Af ~ 2. Perhaps 
the most promising new method is to use chiral perturbation theory. Calculations with only 



the one-loop contribution produce = ^OPn ~ 13.6. Including the next order as 

well as the effects of the A resonance gives [p5| , p^] dp = 10.5 ± 2.0, a„ = 13.4 ± 1.5, 
f3p = 3.5 ± 3.6, f3n = 7.8 ± 3.6. If we compare these two calculations with the Baldin sum 
rule, we find agreement for the predictions for the proton. However, the earlier calculation 
gives better agreement for the neutron than the higher-order one. The most important con- 



tributions to d come from polarizing the pion cloud - about 50-70% p3| . The small value of 
j3 ^ 2 can be attributed to cancellations between the large paramagnetic contribution from 
the A resonance (/?) and the diamagnetic contribution from the pion cloud (A/5) p7| . 

Several experiments have been performed over the past few years with the aim of mea- 
suring nucleon polarizabilities. Most of these have specifically targeted the proton, with a 
differential cross-section for Compton scattering given (to order uj"^) as 

% = % - — \\{dp + 4)(i + cosef + hdp - ^,)(i - cosef] , (4) 

d\l ail Born nip \U! J 12 2 J 



where a is the fine structure constant, and the Born term is that of Powell Note that 
the polarizability terms in equation (^ arise from interference between the polarizability 
amplitude of equation (0) and the Thomson amplitude. This formula shows the sensitivity 
of the cross-section to dp and (3p at different angles. The cross-section at forward angles is 
most sensitive to dp + j3p, while at backward angles only dp — f3p can be measured. At 90° 
the I3p terms drop out completely. The expansion to order uj'^ is only valid to about 100 
MeV [^]. Since the polarizability terms are of this order, too low an energy means that 
the cross-section is does not depend on dp and (3p. Thus the optimal energy range for an 
experiment would be 70-100 MeV, which balances sensitivity with expansion validity. 

A summary of a few proton Compton scattering experiments is given in Table || below. 



The data are those of Moscow experiment [0, Ihinois group ||32|, and Mainz experiment 
More extensive reviews of these and other proton Compton scattering experiments 



can be found in the hterature |T^,|29|. Despite some shortcomings, measurements for the 
polarizabihties in these experiments are in reasonable agreement with each other and with 
theoretical predictions. This is not the case for the neutron. 

A direct neutron Compton scattering experiment would require the use of a neutron 
target to detect a very small cross section. Therefore Coulomb scattering and quasi-elastic 
Compton scattering from a neutron bound in the deuteron have been used. The large 
Coulomb field of a heavy nucleus can induce a large dipole moment. Two separate 1988 
experiments, both using ^°^Pb as the target nucleus, quote large uncertainties: a„ = 12.0 ± 
= 8.0 ± 10.0 A more recent experiment by Schmiedmayer et al. ^6[, also 



10.0, 34 



using a lead target, produced a„ = 12.0 ± 1.5 ± 2.0. This is the value normally quoted for 
the neutron electric polarizability. However, the accuracy of this result has been questioned 
37| , |38[] . It has been claimed that the uncertainties were underestimated, and that the best 



estimate for «„ is really ~ 7 — 19 P^. An experiment to measure quasi- free Compton 
scattering by the neutron bound in the deuteron was carried out by Rose et al. in 
1990, and obtained a„ = ll.Tti^^y.The uncertainties are large, but recent arguments indicate 



that this method is capable of producing more accurate results ^1| and should be revisited. 

It is clear that better experimental measurements of the neutron polarizability are 
needed. This is why we have chosen to investigate deuteron Compton scattering. This idea 



seems to have been originally suggested by Baldin [T^, who calculated the cross-section 
in the impulse approximation. A more extensive calculation was undertaken by Weyrauch 
4^ , |43[| , but certain deficiencies to be discussed below, as well as a lack of emphasis on the 



neutron polarizability, suggest that this question can be reexamined. 
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II. DEUTERON COMPTON SCATTERING CALCULATIONS 



The Feynman diagrams for several of the processes contributing to deuteron Compton 
scattering are shown in Fig. 0. The ovals at either end represent the deuteron wavef unctions, 
while the solid lines are the individual nucleons with which the photons (wavy lines) interact. 
Not all possible combinations are shown; interactions can occur on either nucleon with either 
an incoming or outgoing photon. Figure |I](a) is known as the seagull diagram, arising from 
the term in the Hamiltonian that is proportional to A^. All other one-body interactions, 
which are at least of order u^, are represented by Fig. |T](b). These include the terms which 
depend on the polarizabilities a and (3, as well as some relativistic corrections. Dispersive 
diagrams without meson exchange are depicted in Fig. [l|(c)-(d). At lowest order, each 
vertex can be either an electric or magnetic multipole interaction. Our notation is that 
the intermediate state includes the effects of the un-denoted n-p interaction. Two of the 
meson-exchange effects which are included in this calculation are shown in Figure |l|(e)-(f). 
The term of Figure |l|(f) is called the "vertex correction". 



A. Overview of Calculation 

We begin with an overview of the methods used in the calculation. The differential 
cross-section will be calculated non-relativistically, to second order in the electric charge, 
using Fermi's golden rule. The transition matrix Tfi is then given by 

T,.-{dj,P,,^,\H I P. 1.) + ^ ^./2^^ ^ _ _ p./2^^ ^ + ■ ■ ■ , (5) 

where the ■ ■ ■ represents the crossing term, see Eq. ([Al| ). The notation dn^f) is meant to 
represent all quantum numbers needed to describe the initial (final) deuteron state, except 
for the center-of-mass momentum P, which will be absorbed into a momentum-conserving 
delta function. C represents all possible intermediate states. 

The most important terms we include are represented by a non-relativistic Hamiltonian: 
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where the vector potential is given by 



fc,A=±l 

The operators a^^)^ destroy (create) a photon with momentum k and polarization A. Similarly, 
the pion field operator is given by 

0±(f,) = E ^^/^ («Te^^^-^^ + a^e-^^^-^^) . (8) 

Here, a± destroys a 7r±. The tildes in the Hamiltonian indicate vectors in isospin space. The 
approximate Hamiltonian (0) includes the effects of electric and magnetic dipole interactions. 
As discussed in Sect. IIC below and in Ref. [|lT|, higher multipole moments are included also. 

The transition matrix is proportional to a momentum-conserving delta function. We 
remove this delta function and define the scattering amplitude M-ff 



1 271% 5{Pf + kf- Pi-k 
V y/uJiuJf V 



The differential cross-section is given by 



where Ef = Ui + rrid — is the deuteron mass. The spin-averaged cross-section is 

computed by summing over the final and averaging over the initial polarizations. 



B. Seagull Terms 

Using the interaction Hamiltonian H^'^ = i^A^^Xp) accounts for the seagull diagram. 
The resulting contribution to the transition matrix is given by 
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(11) 



( I ,\ 



M,j 



( 1 



y-M/ M/ - Mi Mi 



) 



t22 

02 I t20 -'2 • 



where r = Xp — Xn, and 



f°° qr 
/ drui>{r)j^{—)ui{r). 
Jo 2 



(12) 



The radial deuteron wavefunction is represented by M02 in which 0,2 refer to the orbital 



angular momentum. We use the "B" wavefunctions of Ref. [0], but the numerical results 
do not depend on this choice. 

The terms that depend on the nucleon polarizabilities are very similar to the seagull 
term. Their transition matrices are given by the following expressions: 



M 



fi 



Y 1 01^ 

J 



(13) 



\-Mf Mfj 



( 1 



y-Mf Mf - Mi Mi 



j22 



-ujfUi\/l27r (ki X ex,) ■ [kf x e^^. 



(-l)^-^^^/(/3, + /3„)x 



1 



1^ 



\-Mf Mfj 



y-Mf Mf - Mi M, 



J 



j22 

2,Mj-M,('?)(-^2 +-^2 ~ 



(14) 
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C. Dispersive Terms 



Most of the remaining terms to be calculated are the dispersive terms, which are second- 
order in the interaction Hamiltonian. The contribution of the dispersive term to the transi- 
tion matrix has the form 

^Msr, ( df, Pf, 7/ I H^""' \C,Pc){C,Pc\ H"^' \ d,, P„ 7. ) . . 

hu;, + E,^-Ec-PS/2m, + ze + ^ ^ 

( df, Pf, 7/ I I C, Pc, 7» lf){C, Pc, li, 7/ I H'-' I d„ Pi, 7, ) 
-huf + Ed, -Ec- Pcl'^md + ie 

in which = — J J(^) ■ ^4(^)0?^^. We use ^ as a dummy variable, while x is a nucleon 
variable and r* is a deuteron variable. The term A is handled using a multipole expansion: 



L=l M=-L 



M2L + 1) 

\ L{L + 1) ^^^^ 



where d\^\{0, —i), —ip) is the Wigner-d function (we use the convention of [^]) , which is 
the overlap of the state | i^A ), rotated by the Euler angles (0, —{}, — f), with the state | LM ). 
To simplify the notation, we define functions $i and $/ by 



1 /2^^^^^e^A:K = V$.(a + ---, (17) 



1 



2^^1e-%-« = V$,(a + ---. (18) 



The gradient terms above correspond to the first term in the brackets of Eq. (W); see also 
Eqs. ( |A2| ) and ( [A3| ) . In the low-energy limit, ~^ ^ ' ^Xii so e<l> is responsible for the 

electric dipole interaction. 

Since A contains both a $ and a non-$ term, each if™* is the sum of two terms and 
then there are four different possibilities for each term of Eq. ([T5|) . The largest term, which 
includes an electric dipole El interaction at both A^7 vertices, is the one for which both 
matrix elements contain $. It is calculated in detail in Appendix A, and certain aspects 
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of it will also be discussed in this section. All terms in which one of the matrix elements 



contain $ are calculated and tabulated in Appendix G of Karakowski's thesis |jTT|. These 
include all magnetic and spin-induced interactions occurring at exactly one vertex. Most of 
the terms that do not contain $ are very small; see Appendix G of Ref. [pH]] . 

The largest second-order terms are the ones in which - / J(0 ■ V$(Oc^^^ is substituted 
for all four if™*'s in Eq. ([TsD . To avoid using an explicit expression for J, we integrate by 
parts and use current conservation: 

V ■ m = ~ [h, piO] , piO = E - ^,)- (19) 

The integral over ^ can then be performed to give 

H'^' = -t[H,{e/h)<!>{x,)], (20) 

where H is now the full Hamiltonian. If we examine only the piece of this commutator 
containing the proton kinetic energy, we see that 



^int,p 



2 



-—Pp-A{xp), (21) 



so the $ commutators contain the operators responsible for El transitions. 

We switch to the center-of-mass variables p = {pp — Pn)/'^ and P = pp + pn when we 
evaluate the commutator with the full Hamiltonian. Defining the "internal" Hamiltonian 

2 

H"^P = + V, we get for the term 

i/int _ __^p . ^^^^ik^■r/2^ikrR _ ■ ^jjnp^ (e/^) (f/2)] 6^^^'^. (22) 



rUp 



The term with the P operator will be called the "recoil correction"; its matrix elements 
are calculated in Appendix B. The substitution Xp = f/2 + R generates an exponential of 
the form e*^'^, which combines with other similar exponentials to create the momentum- 
conserving delta function of Eq. (^). The net result of ignoring the recoil operators, therefore, 
is to replace Xp by f/2 and switch from Tj^ to M.fi- 

The interesting physics is contained in the second term of Eq. (P^. Plugging this into 
Eq. (p!5|), we get 
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M 



fi 



E 



df\ [H^P,^f] \C){C\ [H^P,^^] \di) 



^ ^ huji + Ed^ -Ec- Pc/'^md + ie 
df I [i^"^<l'/] \C){C\ [H^P,^^] \d^)\ 



(23) 



hui + Ed^ -Ec~ Plllmd ]' 

The dimensionless operators = {e/h)(^i{f/2) and $/ = (e/h)^ f{f/2) have been intro- 
duced, and this term has been labeled with the superscript "a" . 

Since | dij ) and | C ) are eigenstates of if"^, the commutators can be expanded and some 
cancellations with the denominator can be made. This generates four terms, after adding 
the crossed term of Eq. ( pTSf ) : 



(24) 



where 



2'm„ 



hujf + 



1 2 



hujj 



2md 



c nuj. 



{df\^f\C){C\^.\d.) 



2md 



+ Ed, -Ec + ie 



E 

c 



{df\^,\C){C\^f\d,) 



-huj 



(huf) 



hujf + ^ — — — huoi + 

2nid 

'-{dj\[[H-^M^f\ + 



f 2md 



1^ + Ed^ - Ec + le' 



2m^ 



( df I I di ), 



(25) 
(26) 

(27) 
(28) 



The first two terms correspond to the Feynman diagram Figure |I](b) and (c), and includes 
the case where both vertices are El interactions. The third term is a small recoil correction, 
calculated in Appendix B. The final term is responsible for compliance to the low-energy 
theorems which result from demanding gauge invariance; see Section |11E| . 

The evaluation of Eqs. (|25| , ^6D while more complex than the one-body terms, is straight- 
forward except for the additional complications of the intermediate states and energy denom- 
inators. We can obtain a compact expression by first examining the algebraic representation 
of the $ operator in the appendix, equations ( [A2| ) and ( |A3|) , and noting that it contains 
distinct angular and radial pieces, which we call O and J, respectively (suppressing the sums 
and related indices here for simplicity). Then we find the J^J] term can be rewritten as 
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■^/' = E E E ( ^'11^/ 1 Of I c){c I I /iiM, ) 

i=0,2/'=0,2 



X 



^0 



rx]^ {r)Ji{r)ui{r)dr, 



(29) 



where the angular quantum numbers (LcScJcMc), are denoted collectively by C. The 
energy denominators have been incorporated in the function Xj'"(r). All of the second-order 
terms can be reduced to the form of Eq. (PU]). We determine X/*"!^) by solving 

d'^ m.„ , Lc{Lc + 1) 



r 



(30) 



where i?( 







, -Efir For V(j{r), we use the Reid93 np potential ^6|. This is 
an ordinary, inhomogeneous, second-order differential equation. The boundary condition is 
that must be an outgoing spherical wave at large distances. 



D. Relativistic Corrections 

The major relativistic correction here is the spin-orbit effect. Corrections arising from 
boosting the final deuteron wavefunction into a moving frame are small here and neglected. 
The origin of the relativistic spin-orbit effect can be understood by examining the magnetic 
term in the Hamiltonian of Eq. (|^): 

H" = Y. "-^^S, ■ B(l,). (31) 

j=n,p j 

The relativistic correction to the magnetic field for an object moving with velocity v is 
B ^ B — V X E, which leads to a relativistic term in the Hamiltonian: 

H-- = -Y. '-^^^^S, ■ {v, X 4) . (32) 

j=n,p 3 

Only the proton gives a correction because of the factor Cj. The effects of Thomas precession 
are included. Since v ^ p/m and E = W{r) ~ r for a central potential, the dot product 
above goes as a ■ L, where L is the orbital angular momentum. If we make the minimal 
substitution p^p — eA, and use E = H^^ becomes 
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e(l + 2Ki) 



tr-RC _ t^^x -r ^rbj 



(33) 



and the resulting correction to the scattering amphtude is given by 

Mff = { dj I ■ (6l^ X exjihcof + hco.)i2K, + l)e-'^^-^^/' \ d, ). (34) 



E. Checking the Calculation 

There are checks which can be performed to help ensure that certain parts of the cal- 
culation are correct. One involves the use of the optical theorem to calculate the total 
cross-section citotfor deuteron photodisintegration, 7(i — > np. With our normalization con- 
ventions, the optical theorem is given by 

An 

cTtot = —\mM^{e = 0), (35) 

in which M.fi is our Compton amplitude. Only the diagram of Fig. lb enters here. How- 
ever, all combinations of interactions at the two vertices must be included. In a 1964 paper, 
Partovi calculated the contributions from each of these interactions at 3 different en- 
ergies. The dominant term is the one in which there is an El interaction at each vertex. 
He called this cross-section "approximation A", and added the other terms one by one in 
successive approximations, until every term was included ("approximation I"). These results 
are reproduced in Table along with the corresponding cross- sect ions extracted from this 
calculation. There is reasonable agreement at all energies; any differences can be attributed 
to changes in the short-range parts of the potentials. The meanings of the approximations 
in Table |I[ Approximation A: El only; Approximation B: A + singlet Ml; Approximation 
C: B + E2; Approximation D: C + triplet Ml; Approximation E: D + spin induced triplet 
Ml; Approximation F: E + spin induced triplet M2; Approximation G: F + spin induced 
triplet El; Approximation H: G + retardation corrections to El; Approximation I: H + all 
other terms. 
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We also compare these photodisintegration cross-sections with experimental data up to 
100 MeV. This is shown in Fig. |^. A vertex correction (see Fig. If) that contributes to this 
process, and is responsible for dominant meson exchange correction [^] is also included in 
a second graph of the range 50-100 MeV. There is good agreement with the data. Our 
photodisintegration cross-section also reproduces the correct threshold behavior, including 
pion-exchange current effect, as described by Brown & Jackson 



Another important check is to ensure that the calculation is gauge invariant. Since the 
relativistic two-photon scattering amplitudes must have the form efT^j^e^, the condition 
imposed by gauge invariance is that k'^Tfj^^e'j = e^T^i.k'j = 0. An important consequence of 



this is what is usually referred to as the low-energy theorem. Friar [^] showed, assuming 
only gauge invariance, that the full Compton amplitude must go as 

^2 



M 



(36) 



as the photon energy goes to zero. This result is non-trivial because it involves the tar- 
get (deuteron) mass; at low-energy the seagull amplitude depends on the proton mass : 

^2 



M 



sg 



rrir. 



[ej ■ ei 



(37) 



Therefore, all other terms, including meson-exchange terms, must exactly cancel half of the 
seagull amplitude at threshold. That this occurs is easy to show, if we ignore meson-exchange 
terms. The only term which survives in the low-energy limit, from Eq. (pHj), is 



M 



aA 



p 



+ 



. p 



m. 



di ), 



(38) 



where we have only included the kinetic energy part of H"^^; the potential energy term 
is cancelled by meson-exchange currents and will be discussed shortly. The amplitude of 
Eq. (|37|) can easily be shown to satisfy Friar's low-energy theorem. Since A ^ e a.t threshold 
(the dipole approximation), the definition of $4 Eqs. ( p^Tf ) and (p!8|) implies that '^'i ^ f| • 
ii, "^z ^ f i ■ Evaluating the commutators of Eq. (|38|) then yields 



M 



a4 



(39) 
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which cancels half of the seagull term as required. Our numerical calculation reproduces this 
result. The differential cross-section in this figure includes all terms except meson-exchange 
terms, and is seen to approach the correct limiting value. This plot shows only the forward 
cross-section, but the low energy theorem is satisfied at all angles. 

We now examine the effects of the meson-exchange currents that were previously ne- 
glected. Since the potential energy portion of the double commutator term (Eq. |38|) has not 
yet been accounted for, we suspect that this must cancel the explicit meson-exchange terms 
at low energy. This is indeed the shown by Arenhovel [^] . Fig. ^ shows the four 

pion-exchange diagrams which contribute and cancel the double-commutator \/Opep ^^qj-t^ 
threshold. These four terms and the double-commutator V^op^p term are described in detail 



in in App. I of [0. Fig. ^ shows a graph of the differential cross-section at 0° which includes 
only these 5 terms. It is exactly zero at threshold. This is true at all angles, not just in the 
forward direction. Gauge invariance has been satisfied without explicitly demanding it in 
the formulation of the problem. The cross-section can be seen to have an dependence, 
which means that the amplitude goes like u"^. 



III. RESULTS AND DISCUSSION 

We now compute the deuteron Compton scattering cross-sections at the three experi- 
mentally relevant three energies of 49 MeV, 69 MeV and 95 MeV. 



A. Effects of Various Terms in Cross-Section 

The first step is to study the contributions to the cross-section from the major interac- 
tions. They will be labeled as follows: SG: seagull term, EM: All interactions, except for 
recoil corrections, arising from the multipole decomposition of the vector potential as de- 
scribed in App. A here and App. G of |jTT|]. These include all electric and magnetic dipoles, 
as well as the double commutator term needed to satisfy the low-energy theorem, vr: All pion 
interactions pictured in Fig |^ which are needed to satisfy the low-energy theorem. These are 
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detailed in App. J of [|TT[. a + P: polarizability terms. Unless otherwise specified, the values 
for the polarizabilities are chosen as = 10.9, a„ = 12.0, Pp = 3.3, /3„ = 2.0, as determined 
in pT|j36|] . These major terms MT are displayed in Fig. 5 and on the top of Fig. 6. The data 
from |I|] is displayed on all graphs at 49 and 69 MeV; no data have yet been published for 95 
MeV. The seagull term alone provides a reasonable description of the data at 49 MeV, and 
does no worse than all the terms together at 69 MeV. At low energy, the seagull term has 
a (1 + cos^ 9) dependence. As the energy becomes higher, this dependence is maintained at 
the forward angles, while the cross-section get increasingly smaller at the backward angles. 
It is the dominant term at all energies under consideration. The EM terms raise the seagull 
cross-section and become increasingly important as the photon energy rises. The vr terms, 
on the other hand, do not have a very large effect, even at 95 MeV. These terms exactly 
cancel at threshold and continue to nearly cancel at higher energies. 

These terms are collectively called the Major Terms (MT). They describe the data rea- 
sonably well. The smaller corrections are: RC: The Relativistic Correction. CMC: Recoil 
Corrections (from the Center-of-Mass operator), as calculated in Appendix B. This also in- 
cludes only the largest correction. VC: Vertex Corrections, such as the one pictured in 
Fig. |I|(f). These are detailed in [|TT|]. These terms produce significant effects only at 95 MeV. 
A detailed examination of the EM shows that multipoles higher than L = 1 are unimportant 
below 100 MeV. 



B. Determining Polarizabilities from Deuteron Cross- Sect ions 

Having displayed the major contributions, we now investigate the dependence of the 
cross-section on the size of the neutron polarizabilities. Figs. |^ ^ show the effects of varying 
the magnitudes of a„ and The a„ term has the same (l+cos^ 9) dependence as the seagull 
term at low energies, while the /5„ term goes like cos 6*, so that the magnetic polarizability 
term gives no contribution at 90°. The largest effects for both terms occur at extreme forward 
and backward angles, while the experimental data points are in the middle. Figs. |^ and 8 
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show the changes to the cross-section over the range of a„ = 12.0 ± 4.0 and /3„ = 2.0 ± 4.0, 
which is a range of variation that is shghtly larger than the errors quoted in ||36|. These 
illustrate the sensitivity to the size error. At 49 MeV the effects of varying the polarizabilities 
between the values of an = 8.0 — 12.0 and /5„ = 2.0 — 6.0, are very small compared with the 
experimental error bars. There are only two data points at 69 MeV, and one of these lies 
above our calculated cross-section. The data point at the greatest angle is also difficult to 
describe at 49 MeV. Other calculations show similar inconsistencies with these two points 
@J|J^. Looking at (Fig. ||), we see that using a lower and a higher /3„ than the accepted 
values provide the best description of these points. 

We must keep in mind that the proton polarizabilities also have uncertainties, and they 
contribute to the deuteron cross-section just as much as the neutron. Another set of graphs 



as described above is given for Up and jSp (Figs. |T0| and 11) , using only the ranges as the 



actual experimental error given in |3T|]. The same trends stand out here, most notably that 
a smaller value of a and a larger value of [3 would provide a better description of the data 
points at the back angles. Using the uncertainties in neutron and proton polarizabilities 
would allow a perfect reproduction of the data. 

The chance to determine polarizabilities looks more promising at 95 MeV, as shown in 
Figs. 8 and 11. The effects of the polarizabilities are larger, and the range of ±4.0 probably 
will extend beyond the experimental error bars at angles such as 50° and 140°. 



C. Sources of Error 

We would like to study the effects of some of the smaller terms more closely. For each 
term we estimate the absolute error in the value a„ (we fix fin = 2.0 in this section) caused by 
neglecting various effects. These are listed in Table |IT1| and shown in some figures. First, we 
examine the gauge invariant set of pion terms. The results in Table |I| show that the value 
a„ would be slightly underestimated at the forward angles and overestimated at backward 
angles, and would not be affected around 120°. The magnitude of this error is about the 
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same as the experimental error in the polarizabihty, so these effects need to be included to 
make an accurate determination of a„. The effect is same at all energies because both the vr 
terms and polarizabihty terms have an lo^ dependence. The recoil (CMC) corrections would 
be difficult to see on any graph, and this is also shown in the table. The largest errors are 
actually at the lowest energies, indicating that these terms are of order lo^ . The small values 
for Aa„ at 95 MeV indicate that these terms can be neglected here. 



The relativistic corrections (RC) (Fig. |T^) are not a "small" correction: even at 49 MeV 
they would have to be included in a calculation to make an accurate determination of 
The amount of error that would be introduced by neglecting these terms is comparable to 
an itself at some angles. 

We turn to the vertex corrections (VC) of Fig. If which are also necessary for an accurate 
estimate of the polarizabihty. Their effect is greater than the effects of the other pion terms 
combined. The errors decrease as energy increases, but even at 95 MeV the errors, particu- 
larly at middle angles, is at least as large as previous errors in polarizabihty determinations. 
At backward angles the effect seems to be independent of energy. See Fig. |TB|. Finally, we 
note that the effect of the uncertainly in the vrA^ coupling constant / on computed values 
of the deuteron Compton scattering are very small. 



D. Tensor-Polarized Deuteron Compton Scattering 



Another interesting observable is the tensor-polarized deuteron Compton scattering 
cross-section. Chen H introduced the idea that the quantity 



^ dcr , , ^ da , , ^ _ , da , , ^ 

= ") - Isi^"' = ^' - 



(40) 



be studied because it isolates the importance of pion exchange terms in effective field theory 
treatments of deuteron Compton scattering. Fig. ^ shows that (^)y, as obtained in our 
complete calculation has only a slight dependence on the polarizabilities, which are our 
primary interest. We display our results to facilitate future comparisons between potential- 
model and effective field theory calculations. In our approach, the contributions of magnetic 
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interactions (MI) and relativistic corrections (RC) are the dominant contributions to 
We define the pion-exchange terms to include any pions that give rise to the tensor force 
(and the / = 2 state) in the deuteron. However, to be more consistent with the calculation 
of we now will include only one pion exchange in the pion terms and no pion exchanges in 
the magnetic terms. Therefore, the Z^-state deuteron wavefunction U2{r) will be set to zero 
in all terms except for the Thomson term and the D-state contributions to this term will be 
grouped with the pion-exchange terms. Defined in this way, the pion-exchange effects turn 
out to be less than about 0.5 nb/sr at all angles. The results for the magnetic interactions 
(MI) are shown in Fig. These interactions include all of the dispersive terms of Eq. (|A1|) 
as well as the seagull term of Eq. (|12D computed without the deuteron D-state. In addition, 
we show the effect of including the deuteron D-state in the calculation of the magnetic terms 
(but not the seagull term). This is a non-negligible effect. Finally, the relativistic correction 
(RC) is added (including the D-state). Since this is a spin-orbit effect, we group this with 
the magnetic interactions. Both the relativistic term and the deuteron D-states have a large 
effect in the calculation; this observable is very sensitive to these individual effects. 



IV. CONCLUSION 

We have presented a calculation of deuteron Compton scattering valid for energies less 
than 100 MeV. All terms that we believe to be important are included. The seagull, lowest- 
order multipoles, and polarizability effects dominate the cross-section. Relativistic correc- 
tions must also be incorporated in the calculation. Although smaller than the preceding 
terms, pion-exchange and vertex corrections are also needed to accurately determine the 
polarizabilities. We conclude that recoil effects and multipoles of second-order and higher 
are negligible in this energy range. We are able to achieve a reasonable agreement with the 
existing data; see Figs. 5,7 and 10. 

This is the most extensive study so far of the effects of various terms in the cross- 
section on the value of a„; see Tables III and IV. However, certain corrections which are 
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believed to be small have been omitted. No contributions from the A resonance have been 
explicitly included, although the magnetic polarizability itself is thought to contain a large 
paramagnetic contribution from the A. The A resonance has been shown to have a small 



effect on deuteron photodisintegration at energies less than 100 MeV In addition, 

considering the rather small contributions from vr-exchange currents, heavier mesons such as 
the p have also been neglected. Higher-order relativistic corrections have not been included, 
but may need to be investigated further at energies near 100 MeV because of the large effects 
of the leading-order terms there. 

The only fully detailed calculation of this type that has been published to date was by 
Weyrauch |^ , |43[| . There are several important differences between his calculation and the 
present one. One of these is the A^A^ potential used. He used a separable A^A^ T matrix, while 
here we have found the Green's function for an intermediate A^A^ state interacting via the 
Reid93 potential. He also assumes meson-exchange currents and relativistic corrections to be 
negligible. Probably the most important difference, however, is in gauge invariance. Gauge 
invariance is broken in the separable potential model and needs to be artificially restored. 
We have made no such assumptions here. Gauge invariance is a natural consequence of 
including all of the appropriate terms. 

The cross-section predicted by Weyrauch is higher than the experimental data. How- 
ever, recent results of Arenhovel 0] and Levchuk and L'vov are in qualitative agreement 
with the data. Their work is similar to ours in that they also take diagrammatic approach 
using realistic potentials and including meson-exchange currents. The full details of their 
calculations are not yet available, so precise comparisons are not easily made. Lastly, a new 
calculation using effective field theories is also in qualitative agreement with the data [§. 
We have also followed |^ in studying tensor-polarized Compton scattering. In our formu- 
lation the effects of the deuteron D-state and the relativistic correction are the dominant 
contribution. 

We return to our basic question of the feasibility of making an accurate determination of 
the neutron polarizability from a deuteron Compton scattering experiment. The available 
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data at 49 and 69 MeV show that the ranges = 12.0±4.0 and (3n — 2.0±4.0 are reasonable, 
but a more accurate estimate is difficult since the polarizability terms are small relative to 
the experimental error bars. The uncertainties in the proton polarizabilities must also be 
taken into account. These problems are compounded by the data point at the greatest angle 
at each energy, which we all authors find difficult to describe. The closest description of all 
of the experimental data comes from lowering the value of q;„ and raising 

Perhaps the most promising estimates of the polarizability will come from data at higher 
energies. Assuming that the experimental error bars are not significantly larger, we should 
find that the range of polarizabilities that fit the data is smaller. Sensitivity to both a and 
P is largest at the forward and backward angles, while the cross-section at 90° is the least 
sensitive to the polarizabilities. Many of the "small" corrections, however, are particularly 
small at 120 — 140°, and the cross-section is still very sensitive to the polarizabilities at these 
angles. The proton polarizabilities, unfortunately, will also make a larger contribution to the 
cross-section, and there is no way to determine the neutron polarizabilities more accurately 
than those of the proton in this experiment. This energy is at the limit of the validity of 
the low-energy expansion, which may introduce errors as well. In any case, studies at higher 
energies should provide additional information on the polarizabilities. 
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TABLES 



TABLE I. Experimental results for Up and f3p 



Moscow(1975) 



a;i(MeV) 



70-110 



^(degrees) 
90,150 



result 
ap + j3p = 5.8 ± 3.3 
ap-Pp = 17.8 ± 2.0 



Illinois (1991) 



32-72 



60,135 



ap = 10.9 ±2.2 ±1.3 



/3p = 3.3 =F 2.2 =F 1.3 



Mainz(1992) 



i, 132 



180 



Pp — ' -^^-2. 37-2. 1 



05 



ap = 10.62_i,i9_i,o3 



/3p = 3.58+i 



+1.25+1.07 
19-1.03 



Illinois/SAL(1995) 



70-148 



90,135 



ap = 12.5 ±0.6 ±0.7 ±0.5 



= 1.7 T 0.6 T 0.7 T 0.5 
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TABLE II. Comparison of photodisintegration cross-sections calculated in several approxima- 
tions, ap is the cross-section of Partovi |^^, while ax is from this calculation. 



Approximation 20 MeV 80 MeV 140 MeV 





ap {fib/sr) 


ax ifJ-b/sr) 


ap (fib/sr) 


ax ifib/sr) 


ap (fib/sr) 


ax {fJ-b/sr 


A 


579.1 


583.3 


77.15 


80.54 


34.04 


34.56 


B 


589.2 


593.2 


83.50 


85.83 


39.77 


38.59 


C 


591.5 


595.0 


84.55 


86.86 


40.38 


38.99 


D 


591.7 


595.1 


84.85 


87.06 


40.67 


39.16 


E 


591.6 


594.7 


84.73 


86.61 


40.55 


38.70 


F 


592.1 


595.3 


85.43 


87.68 


41.31 


40.26 


G 


591.9 


594.6 


90.45 


90.22 


47.07 


43.33 


H 


588.2 


591.2 


87.52 


86.36 


44.64 


39.44 


I 


588.2 


591.2 


87.40 


86.40 


44.53 


39.52 
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TABLE III. Error in electric polarizability generated by omitting the gauge-invariant set of 
pion terms tt and the recoil corrections (CMC) 



49 MeV 

-2.9 
-2.7 
-2.1 
-1.2 

0.0 

1.1 

1.5 



TT 

69 MeV 

-2.6 
-2.4 
-2.0 
-1.4 
-0.3 

0.9 

1.4 



95 MeV 

-2.0 
-1.9 
-1.6 
-1.7 
-0.9 

0.5 

1.1 



49 MeV 

0.0 
-0.5 
-1.1 
0.1 
1.4 
0.6 
0.0 



CMC 
69 MeV 

0.0 
-0.2 
-0.5 

0.2 

0.8 

0.3 

0.0 



95 MeV 

0.0 
-0.1 
-0.2 

0.3 

0.5 

0.2 

0.0 



TABLE IV. Error in electric polarizability generated by omitting the relativistic correction 
(RC) and vertex correction (VC) 







RC 






VC 




^(deg) 


49 MeV 


69 MeV 


95 MeV 


49 MeV 


69 MeV 


95 MeV 





-7.9 


-9.2 


-11.9 


3.7 


1.4 


0.9 


30 


-8.7 


-10.2 


-13.2 


3.9 


1.9 


1.7 


60 


-9.7 


-11.8 


-15.8 


3.7 


2.7 


4.1 


90 


-4.9 


-7.1 


-10.1 


-1.1 


1.2 


5.0 


120 


5.1 


4.1 


3.4 


-8.0 


-3.4 


0.9 


150 


10.0 


9.8 


10.3 


-10.6 


-5.9 


-2.4 


180 


11.1 


11.1 


12.0 


-11.1 


-6.6 


-3.3 
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FIGURES 
FIG. 1. Feynman diagrams 

FIG. 2. Total photodisintegration cross-sections compared with experimental data. 

FIG. 3. Gauge invariant set of pion-exchange diagrams 

FIG. 4. Forward differential cross-section including only pion-exchange and dou- 

ble-commutator contributions 

FIG. 5. Contributions of different terms. 

FIG. 6. Contributions of different terms at 95 MeV. Both major terms and smaller corrections 
are displayed 

FIG. 7. Effect of varying and separately. /3„ = 2.0 in the top graph and = 12.0 in 
the bottom. 

FIG. 8. Effect of varying a„ and separately at 95 MeV. /3„ are as in Fig. 7 
FIG. 9. Effect of varying both q;„ and /3„ . 
FIG. 10. Effect of varying both ap and 

FIG. 11. Sensitivity to values of a,/? at 95 MeV 

FIG. 12. Effect of neglecting relativistic terms RC on the differential cross-section. The full 
calculation is labelled FC. 
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FIG. 13. Effect of neglecting vertex corrections VC on the differential cross-section 

FIG. 14. Tensor-polarized deuteron Compton scattering cross-section, with and without polar- 
izabilities. The values of a„ = 12.0, Pn = 2.0, Op = 10.9, and /Sp = 3.3 are used in the solid curve. 
All other interactions are included. 

FIG. 15. Tensor-polarized deuteron Compton scattering cross-section, including only magnetic 
interactions(MI). The dotted curve contains all of the dispersive terms, but only the contributions 
from the deuteron S'-state are included. These contributions are added in the dashed curve. Rela- 
tivistic corrections (RC) are added in the solid curve. 



APPENDIX A: LEADING-ORDER DISPERSIVE TERMS 

The remaining terms will be calculated using second-order perturbation theory: 

^disp _ V- [ ( df,Pf,7f I H'-' I C,Pc ){C,Pc I H'-' I ci„P„7. ) 

nu, + E,^-Ec-Pl/2m, + ie ^ ^ > 

{ df, P/, 7/ I I C, Pc, 7., lf){C, Pc, 7i, 7/ I H'-' I d., Pj, l^ ) 
-hujf + Ed, -Ec- Pcl'^md + is 

Here C denotes the internal quantum numbers of the intermediate np state, and X^c 

shorthand for all sums and integrals which are needed to describe this complete set of states 

except for the center-of-mass states. These are written separately as Pc- At this point, 

it is most convenient to use if™* = — / J(^) ■ A{S)d?^-, and then expand A into multipoles 

according to equation (|16|). There are 3 terms in this expansion; therefore, equation (|A1|) 

contains a total of 18 terms. The largest contributions at low energies (which we denote as 

M"fi) arise from the gradient operator in Eq. (jlBj). These are the leading-order contributions; 

the other terms are of order uj/rriN smaller. This gradient term includes the case where there 

is an El interaction at both 7A^ vertices. 

With this in mind, we define $j(r) and by 
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1 + j jLii^ir)YLMir) 



1 O-rrfi °° ^ 



27r(2L' + 1) / d\ , , 



which means that 



1 iT^ri. 



After a lengthy calculation [|TT| one finds 



(A3) 



(A4) 
(A5) 



(A6) 



where 



' 2 oo 1+L Jc+Sc 

E E E E E Ex 

- 1=0,2 l'=Q,2 L,L'=1 Sc=0,l Jc=\l-L\ Lc = \Jc-Sc\ 



{Tiuji){fiujf) [ 2md 



(2L + 1)(2L^ + 1) 
\LL'{L + l){L' + 1) 



V 



J< 



c 



X 



Jc L 

-A,; - Mi \i Mi 



( mil Yl, WLcScJc ){LcScJc\\ Yl ||/11 ) x 



dr r M,(r)x/^(r)^/'i(^). 



uJir, 



We have used the definitions Eq = hoji — -^^^^^ and ip^i.^) = Jl^x) + x-^ii^{x). 

' oo 1+L Jc+Sc 

E E E E E Ex 

/=0,2 r=0,2 L,L'=1 5c=0,l Jc=\l-L\ Lc = \Jc-Sc\ 



(A7) 



-27re2 



{huji){hujf) 



hujf + 



2m. 



d 



,L+L'^_^^L'-Xj+Jc-X,-^{Lp^ 



^ 1 L J. \^ 



y-Mf \i -\ + Mf J 



f-Mi-Xi-Xf 

Jc L' 1 

A, - Mf Mf- Mi- \i M, 



(2L + 1)(2L^ + 1) 
^ LU{L + l){U + l) 

\ 



X 
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( Yl WLcScJc ){LcScJc\\ Yu IKH ) X 



(A8) 



where Eq = —huf 



M 



a3 



2md 



2md 



L+L' oo oo 
J 1=0,2 V =0,2 l=\L_Li\L=l L'=l 



J-^Mf-M,-\-\f 



(O,-^,-0)(2L + l)(2L'+l) 



2L + 1 







L V L 

\i Mf - Mi- \i Mi- M 



U 1 L A 



X 



Mf Mf - Mi Mi 



X 



(A9) 



The final term involves the double commutator of the Hamiltonian. The kinetic energy 
contribution is given here and potential energy below. 
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The other dispersive terms are listed in [11 



APPENDIX B: RECOIL CORRECTIONS 



We will now return to the recoil corrections to the transition matrices which arise from the 
P operators in Eqs. (|15| )and (|2^ ) The main contributions from these terms are 
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Since we are working in the lab frame, P | P, ) = 0. Wc integrate the center-of-mass variables, 
use A — V$, and evaluate the commutator so that 7^°°'"'^^'' becomes:. 
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We list the results for the other Alj™ terms: 
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APPENDIX C: VERTEX CORRECTIONS 

We calculate terms, such as the one shown in Fig. If,, in which a (7,7r) reaction on 
one nucleon is followed by the absorption of the pion on the second. This is followed (or 
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preceeded) by a photon emission from the second nucleon. These are vertex corrections to 
the "ordinary" second-order diagrams. 
The result is 
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where the Green's function G^^ {r,r'; Eq) is defined by 



(CI) 



Eo-H^' + ie 



r). 



The other correction terms are: 
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FIGURES 




(f) 7/ 
FIG. 1. Feynman diagrams 




FIG. 2. Total photodisintegration cross-sections compared with experimental data. 
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(a) 7/ (b) 




(c) 7/ (d) 7/ 



FIG. 3. Gauge-invariant set of pion-exchange diagrams 
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FIG. 4. Forward differential cross-section including only pion-exchange and double commutator 
yOPEP contributions 
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FIG. 5. Contributions of different terms. 
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FIG. 6. Contributions of different terms at 95 MeV. Both major terms and smaller corrections 
are displayed 
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FIG. 7. Effect of varying q:„ and /3„ separately. /3„ = 2.0 in the top graph and q:„ = 12.0 in 
the bottom. 
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FIG. 8. Effect of varying and /3„ separately at 95 MeV. an,Pn are as in Fig. 7 
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FIG. 9. Effect of varying both q:„ and /3„ 
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FIG. 10. Effect of varying both ap and /3, 
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FIG. 11. Sensitivity to values of a,/? at 95 MeV 
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FIG. 12. Effect of neglecting relativistic terms on the differential cross-section. Our full calcu- 
laiton is labelled FC 
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FIG. 13. Effect of neglecting vertex corrections VC on the differential cross-section 
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FIG. 14. Tensor-polarized deuteron Compton scattering cross-section, with and without polar- 
izabiUties. The values of Q!„ = 12.0, /3„ = 2.0, ap = 10.9, and (3p = 3.3 are used in the solid curve. 



All other interactions are included. 
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FIG. 15. Tensor-polarized deuteron Compton scattering cross-section, including only magnetic 
interactions(Ml). The dotted curve contains all of the dispersive terms, , but only the contribu- 
tions from the deuteron 5'-state are included. These contributions are added in the dashed curve. 
Relativistic corrections (RC) are added in the solid curve. 
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